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I .  INTRODUCTION 


1.1  Background  .  The  peak  overpressure  associated  with  a  blast  wave, 
propagating  from  an  accidental  explosion  In  an  underground  munition  storage 
stte.  Is  the  damaging  mechanism  that  governs  the  distance  at  which  Inhabited 
buildings  may  be  located.  There  Is  a  range  of  peak,  overpressure,  from  50 
millibars  (0.725  osi)  to  86  millibars  (1.2  psl)  ,  which  has  been  established  as 
the  criterion  for  acceptable  damage  to  an  Inhabited  building.  NATO  countries, 
In  general,  use  the  50  millibars  while  the  United  States  use  86  millibars. 
There  ate  also  different  methods  used  to  predict  the  distance  one  might  expect 
these  peak  overpressures.  These  differences  In  the  peak  overpressure  for 
acceptable  damage  and  the  methods  for  predicting  the  distance  at  which  this 
pressure  would  occur  are  of  primary  Interest  to  this  report. 

1.2  Ob  joctlves .  The  objectives  of  this  study  are  to  determine  the 
rationale  for  current  criteria  for  both  the  U  .8 .  and  NATO  countries,  to  assess 
weaknesses  In  the  different  approaches,  and  to  establish  a  new  recommendation 
based  on  scientific  experiments  and  theoretical  calculations. 

2 .  RESULTS 

2.1  Literature  Search .  An  extensive  literature  search  was  made  and  a 
total  of  24  reports  reviewed  In  detail.  These  are  listed  as  References  1-24. 
These  reports  Included  small  scale  and  shock  tube  experiments,  and  computer 
calculations  . 

2.2  Chamber  Pressure  and  Exit  Pressure  .  The  various  parameters  that 
govern  the  blast  propagation  outside  of  an  underground  tunnel  are  the  storage 
chamber  dimensions  and  volume,  passageway  dimensions  and  volume,  mass  and  type 
of  explosive  stored,  exit  pressure,  tunnel  diameter,  and  the  angle  off  of  the 
zero-degree  axis. 

The  mass  of  explosive  and  volume  of  the  storage  chamber  ate  needed  to 

establish  the  loading  density.  One  of  the  Norwegian  reports'*  concentrated  on 
the  bulld-up  of  pressure  In  the  storage  chamber  by  measuring  the  pressure 
versus  time  for  different  loading  densities,  types  of  explosive,  and  vent 
areas.  The  experimental  results  compared  quite  well  with  the  output  from 

Proctor's  INBLAST  computer  code  Although  the  chamber  pressure  is  one  of  the 
Important  parameters  and  depends  on  loading  density,  early  equations, 
developed  to  predict  the  exit  pressure  from  the  tunnel,  used  loading  density 

rather  than  chamber  pressure  J  The  equation  established  for  predicting  the 
exit  pressure  Is  approximated  by 

Pw  24(Q/Vt)°’66,  (l) 


1 


where 


exit  pressure,  bars 


o 

"  w 


Q  3  exolosive  mass,  kg 

3 

V  =  total  volume,  m  . 


In  English  units,  the  equation  becomes 


where 


P  3  2172  (W/V  )0  .66  , 
w  t 


P  *  exit  pressure  In  psl, 
w 

W  =  explosive  mass  in  lbs, 


(2) 


and 


V  3  total  volume  in  ft  . 


When  the  passageway  or  exit  tunnel  cross-section  Is  smaller  than  the 
chamber  cross-section,  then  an  attenuation  of  the  shock  was  considered  and 

other  equations  for  P  were  developed.7 


w 


where 


and 


P  3  17.1  (Q/V  )0'607  (A  M  )0*19’ 
”  t  j  c 


w 


3  area  of  exit  tunnel, 

A^  3  area  of  storage  chamber  at  exit 


w 


(See  Figure  l), 
bar,  Q  *  kg,  and  V 


In  English  units,  the  equation  becomes 


(33 


where 


P  3  943  (W/V  r607  (A, /A  )°'19, 
w  t  3  c 

3 

P  3  psl,  W  *  lbs,  and  V  *  ft  . 
w  t 


(4) 


In  Reference  17,  a  new  equation  was  developed  to  predict  the  exit  pressure 

P  3  16.4  (Q/V  3°*54  (A, /A  )0*24*  (5) 

w  t  3  c 


P  .  This  equation  Is: 


As  can  be  seen,  this  Is  a  variation  of  Equation  3.  A  comparison  of  Equations 

3  and  5  shows  that  at  the  lower  loading  densities.  Equation  5  predicts  higher 

values  for  P  ,  while  at  higher  loading  densities  (Q/V  >  303,  Equation  3 
w  t 

predicts  higher  values  of  ?  . 
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A  method  developed 

at  URL  considers  the 

total  volume  pressure 

'vt 
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governing  parameter  rather  than  the 
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predict  the  chamber  pressure  for  a  qiven  explosive  and  storage  density.  In 
the  BKL  method,  the  total  volume  is  used  in  the  equation  and  the  sane 
attenuation  factor  usinq  A./A^  also  appears  in  the  equation,  as  follows: 


where 


P 

w 


(P..r)'l,H5  (A. /A 
Vt  |  c 


P  and  P,  are  in  bars, 
w  Vt 


(h) 


In  F.nqlish  units,  this  translates  to: 


P 

w 


1.711  (Pvt)°-H1  (A./A(.)0*'9 


(7) 


where  P  and  P.,  arc  in  psi . 

w  Vt 

q 

A  plot  of  H  (psi)  versus  w/v^  (Ih/ft  )  is  presented  in  Fiqure  2  for  both 

TNT  and  PF.TN.  This  is  to  illustrate  that  the  chamber  pressure  for  each 
specific  explosive  should  be  calculated  rather  than  usim*  a  TNT  equivalence 
factor.  In  this  illustration,  PKTN  shows  a  lower  efficiency  than  TNT  at  the 
low  loading  densities,  but  becomes  higher  above  a  loading  density  of  O.OM 

Ibs/tt*.  The  total  volume  pressures  as  a  function  of  loading  density  for 
various  explosives  are  listed  in  Table  I.  This  table  was  taken  from 
Reference  19. 

When  Equation  6  is  compared  with  Equation  ‘i ,  the  values  of  the  predicted 
exit  pressures  for  Equation  ft  are  larger  at  the  higher  loading  densities. 
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Equation  5 

Equation  6 

Q/vt 

P  ,  bar 
w 

P  ,  1 
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This  increase  in  t. he  value  of  the  exit  pressure  may  be  justified  because  it 
can  be  seen  in  Figure  2  that  as  t  lie  loading  density  increases,  the  chamber 
pressure,  P  ,  increases  quite,  rapidly. 

2.3  Outside  Pressure.  A  metliod  for  predicting  the  pressure  propagating 
outside  of  the  tunnel  exit  and  along  different  radials  was  developed  and 
presented  in  Reference  10.  The  basic  equation  is  presented  as  follows: 


where 


and 


AP/P. 

AP 

P 

w 

R 

0 

t 

9 


I  (R/0f  *  ’S/  I  +  (0/56)^  , 

pressure  at  target  in  bar  or  psi  , 
exiL  pressure  in  bar  or  psi, 
distance  to  target  in  m  or  ft, 
tunneL  diameter  in  m  or  ft, 
angle  in  degrees,  off  zero  axis. 


(«) 


Equation  8  has  been  plotted  in  Figure  1,  along  with  data  points  taken  from 
experiments  reported  in  References  10,  11,  and  14-17.  It  is  interesting  to 
note  that  data  from  References  II  and  14-16  were  generated  from  shock  waves 
exiting  from  shock  tubes. 


In  practical  use,  the  desired  parameter  is  the  distance  R  at  which  a 
selected  pressure  would  occur.  Therefore,  Equation  8  may  be  rewritten  as: 
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Figure  3.  A  P/Py  versus  R/Df  Along  Different  Radials 


These  equations  were  developed  by  the  Norwegians  and  presented  In 
Reference  10.  The  distance  Rq  along  the  zero  line  for  a  given  pressure  may  be 

multiplied  by  the  attenuation  factor, 

-  ■  [■  *  wr "  ■ 

to  obtain  the  distance  along  any  radial  at  which  the  same  given  pressure  might 
be  expected.  This  attenuation  factor  AF  Is  plotted  versus  angle  off  the  zero 
axis  in  Figure  4.  The  dashed  lines  in  Figure  4  represent  the  present 
attenuation  system  where  sectors  are  used  rather  than  a  continuous 
attenuation . 


A  second  method  was  proposed  by  the  Norwegians.  In  this  report, 
equation  in  the  form  of  Equation  8  was  presented  . 

AP/P  -  (1  .2987  R/D„)~1,29870cn),  (10) 

w  c 

where  kn  Is  an  attenuation  factor  for  different  sectors.  0°-30°:  kn  -  1  and 
30°-60°:  kn  =»  0.74. 

If  we  put  Equation  10  into  the  form  of  Equation  9,  then  we  have: 

R  -  D,  (0.77)  (P  /AP)°‘77(kn).  (II) 

t  w 

When  values  of  R/D  andAP/P  from  this  equation  are  compared  with  Figure  3, 

t  w 

they  fall  below  the  curve  established  for  Equation  8.  The  attenuation  factors 
for  distance  k  versus  angle  sectors  are  0-30:  kn  «  l,  30-60:  kn  -  0.89, 
60-90:  kn  -  0.67,  90-120:  kn  -  0.5,  and  120-180:  kn  -  0.25. 

Equations  10  and  11  were  developed  from  small  scale  experiments  and  the 
data  falls  along  the  calculated  curve,  but  it  is  recommended  in  this  report 
that  Equation  9  be  used  to  calculate  the  distance  at  which  selected  peak 
overpressures  should  occur  .  This  recommendation  is  based  on  the  fit  of  data 
from  other  sources  as  shown  in  Figure  3. 

2.4  Other  Methods  Considered.  There  are  two  other  methods  that  were 
proposed  for  consideration  as  criteria  for  predicting  the  distance  at  which  an 
inhabited  building  could  be  located . 

20 

The  first  method  was  submitted  by  the  Norwegians.  The  basic  equation 
to  predict  the  distance  to  expect  a  peak  shock  pressure  of  50  mbar  is  as 
follows : 

R  -  18.8  (Q/V)°*265  (Q/nk)0*283,  (12) 

o 
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where 


Q  ■  explosive  mass  In  kilograms, 

3 

V  *  volume  of  storage  chamber,  m  , 
n  »  l  when  storage  site  has  only  one  exit,  or  when 
there  ate  mote  than  one  and  the  blast  waves 


interact . 

n  *  2  when  there  are  more  than  two  exits  and  the 
blast  waves  are  not  expected  to  Interact, 
k  3  3  If  the  branch  passageway  between  the  storage 
chamber  and  the  main  passageway  has  the 
following  characteristics: 

-  crosssectional  area  is  not  greater  than  1/2 
the  main  passageway  area, 

-  length  is  not  less  than  2/3  of  the  requited 
interval,  and 

-  the  angle  between  main  passageway  and  branch 
passageway  Is  within  the  interval  of  60°  to 
120°. 

k  ■  1  for  all  other  cases . 

Equation  12  covets  the  section  0°  to  30°.  For  sector  30°  to  50°,  the 

constant  18.3  is  reduced  to  16.9;  from  60°  to  90°,  18.3  becomes  12.5;  from  90° 

to  120°,  18.8  becomes  3.1;  and  for  120°  to  130°,  the  constant  13.3  is  reduced 
to  4.7.  These  attenuation  factors  for  distance  are  the  same  as  the  dashed 
iine9  in  Figure  4.  A  comparison  of  this  method  with  other  methods  will  be 
presented  later  in  this  report. 
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gure  4.  Attenuation  Factor  versus  Degrees  Off-Axis 


A  second  method  which  Is  quite  similar  to  the  one  just  discussed  was 

oi 

proposed  by  Paul  Price,  DOD  Explosives  Safety  Board.  This  equation  Is 
presented  as  follows: 


where 


and 


R 

o 


F 

y 


_  0 .265  ,  .X 

Fy  (W/V)  (W/nk)  y, 

Is  a  function  of  explosive  mass. 

Is  a  function  of  explosive  mass, 


W  =  explosive  mass  In  pounds,  lbs. 


V  =  volume  of  storage  chamber,  ft  . 


(13) 


n  and  k  have  the  same  definition  as  given  In  the  previous  method. 

When  calculating  R  for  various  charge  masses,  these  factors  are 


listed  below: 

0-100000 

o 

100000-250000 

250000-500000 

F 

X 

F  X 

F  X 

y 

y 

y  y 

y  y 

02 

0  .233 

5  .29  0.531 

115  0 .233 

This  calculation  for  R  includes  the  0°  to  50°  sector.  For  the  other 

o 

sectors,  use  the  dashed  line  attenuation  factors  given  in  Figure  4. 

Both  of  these  methods  have  certain  requirements  which  must  be  met.  The 
first  method  states  that  the  cross-section  of  the  main  passageway  must  not  be 

larger  in  cross-section  than  20  m',  the  tunnel  roughness  must  be  at  least  5%, 
and  the  length  of  the  passageway  must  be  at  least  100  meters.  There  are  no 
corrections  given  for  smaller  area  tunnels,  shorter  tunnels,  or  longer  tunnels. 

The  method  proposed  by  Paul  Price  does  not  specify  tunnel  cross- 

sectional  area,  but  states  that  if  the  tunnel  is  longer  than  330  feet,  then 

reduce  the  distance  R  by  23%. 

o 

A  third  method  for  comparison  is  the  current  one  in  Reference  22. 

Distances  in  this  standard  are  based  on  the  equation: 

Ro  “  76<V1/3’  (l4) 
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where 


W  =  W/nk . 
r 

n  and  k  are  similar  to  previous  description. 

=  range,  ft,  and 

W  =  explosive  stored,  lbs. 

The  loading  density,  chamber  volume,  and  passageway  length  or  diameter  are 
not  required  for  this  method  of  calculation. 

3 .  COMPARISON  OF  RESULTS 

3  . 1  Descr iption  of  Method . 

22 

3.1.1  Method  1 .  This  method  is  published  in  the  current  safety  manual 
and  will  be  presented  to  show  that  in  most  test  cases,  it  is  very 
conservative.  Equation  14  is  used  in  Method  1. 

I  /  3 

R  (ft)  =  7<S  (w  /nk)  '  , 
o  r  ’ 

where  R  is  the  inhabited  building  distance  for  the  0°-3O°  sector  (1.2  psi  or 
o 

S3  mbat )  . 

3.1.2  Method  2 .  Method  2  was  proposed  by  the  Norwegians  in  Reference  20. 
In  Method  2,  Equation  12  is  used  as: 

R  (m)  =*  IS.S  (Q/V)°'265(Q/nk)°'2S3 

o 

This  method  was  detailed  in  Section  2.4  above. 

3.1.3  Method  3 .  This  method  is  similar  to  Method  2,  with  the  exception 
of  a  change  in  constants  and  exponents  depending  on  the  change  in  mass.  In 
Method  3,  Equation  13  is  used. 

R  (ft)  =  92  (W/V)°'265(W/nk)0,283 
o 

This  equation  is  used  for  a  W  of  0  to  100000  lbs.  Here  again,  Rq  applies  to 
the  0°  to  30°  sector  . 

3.1.4  Method  4 .  This  method  is  one  proposed  in  Reference  7.  It  requires 

the  geometry  of  the  storage  site  and  mass  of  explosive  in  order  to  calculate 

the  exit  pressure,  P  ,  and  a  second  equation  is  used  to  calculate  R  .  The 

w’  o 

first,  Equation  4,  is 

P  (psi)  =  943  (W/V  )°,607(A  /A  )°'19, 

W  C  j  c 
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then  Equation  9  in  English  units  is: 


R  (ft.)  =  !).„  (  l  .1  7  D(AP/P  )  , 

U  I  w 

then  attenuation  factors  are  applied  for  the  different  radiais  as  presented  in 
Figure  4. 

3.1.5  Method  5.  This  method  was  also  developed  by  the  Norwegians.*7 
Equation  5  in  English  units  becomes: 

Pw  (psi)  =  1064  (W/VL)°*54(Ai/Ac)0,24; 

then  Equation  11  becomes: 

R  (ft)  =  I)  (0.77)(P  /AP)0‘77 

O  L  W 

where  R  is  used  for  the  0°  to  30°  sector.  For  30°-60°,  use  0.89  R  ,  60°-90° 
o  o 

use  0.67  R  ,  for  90°-120°  use  0.50  R  ,  and  for  120°  to  180°  use  0.25  R  . 

o  o  o 

3.1.6  Method  6.  This  method  was  developed  at  URL,  and  is  being  proposed 
as  a  new  criterion  for  predicting  the  distance  at  which  a  specific  peak 
overpressure  should  occur.  The  major  difference  in  this  method  is  that  the 
pressure  in  the  overall  chamber  and  tunnel  volume  is  used  in  Equation  7  rather 
than  loading  density. 

A  O']  IQ 

Pw  (psi)  =  1.733  (Pyt)  *  ’(A./Ac)0*J  . 


Then  Equation  9,  in  English  units,  is: 

R  (ft)  =  D  (1.173HAP/P  )”°*74. 
o  t  w 

3.2  Comparision  of  Methods.  A  comparision  of  the  six  methods  will  be 

made  where  the  initial  storage  site  parameters  are  the  sane,  so  that  a  direct 

comparison  can  be  made.  We  will  assume  there  is  only  one  tunnel  exit,  then  n 

*  1  and  the  criteria  are  met  to  make  k  *  1 .  The  exit  tunnel  diameter  is  16.6 

feet  and  calculations  will  be  made  for  the  distance  to  1.2  psi  and  0.725  psi 

(50  mbar).  The  ratio  A. /A  from  Figure  1  is  0.23. 

J  c 
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3.2.1  Comparlsion  of  Six  Methods  -  Increase  In  Charge  Mass.  In  Table  2, 
the  volume  of  the  storage  chamber  and  the  passageway  tunnel  remained  constant 
while  the  amount  of  explosive  was  increased  from  2204  lbs  to  11020  lbs,  an 
increase  of  five  times.  The  increase  in  distance  ranged  from  a  factor  of  1.71 
to  2.41.  With  the  exception  of  Method  1,  the  spread  of  distances  for  the  five 
other  methods  is  within  +  11%. 


TABLE  2.  Comparison  of  Six  Methods  -  Increase  in  Charge  Mass 


Loading 


Charge 

Den si ty 

0.725 

Mass  W 

W/V 

C  3 

(  lbs/ ft'*) 

w/v 

(Ibs/ft S 

l .20  psl 

psi 

P 

w 

Method 

(  lbs) 

R  -ft 
— o - 

R  -ft 
— o - 

psi 

1 

2204 

0.062 

0.021 

980 

— 

— 

2 

2204 

0.062 

0.02  1 

— 

4  3b 

— 

3 

2204 

0.062 

0.021 

OS  9 

— 

— 

4 

2204 

0.062 

0.021 

3,36 

561 

63 

5 

2204 

0.062 

0.021 

364 

637  ' 

93 

6 

2204 

0.06  2 

0.021 

444 

644 

32 

1 

11020 

0.312 

0.105 

1691 

— 

— 

2 

11020 

0.312 

0. 105 

— 

1052 

— 

3 

11020 

0.312 

0.105 

941 

— 

— 

4 

11020 

0.312 

0.105 

797 

1157 

181 

5 

11020 

0.312 

0.105 

709 

1045 

221 

6 

11020 

0.312 

0.105 

800 

1 162 

182 

NOTE:  Storage  site  dimensions  constant. 
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3.2.2.  Increase  in  Chamber  Volume  and  Explosive  Mass.  Tn  Table  3,  the 
amount  of  explosive  was  increased  by  n  factor  of  10,  and  the  chamber  volume 
was  increased  by  a  factor  of  10,  so  the  loading  density  remained  the  same 
(0.624).  The  volume  of  the  tunnel  passageway  was  increased  approximately  30%. 
This  changed  the  loading  density  of  the  total  volume  from  0.211  to  0.499.  The 
distances  calculated  for  1.2  and  0.725  psi  at  the  0.211  loading  density  are 
within  +  7%  with  the  exception  of  Method  I.  When  the  loading  density  of  the 
total  volume  was  changed  to  0.499,  the  spread  of  distances  increased  to  +  20%. 
Methods  4-6  are  usually  quite  consistent  in  that  Method  5  calculates  values 
that  are  less  than  the  other  two. 


TABLK  3.  Comparison  of  Six  Methods  -  Increase  in  Charge  Mass 


and 

Charge 
Mass  W 

Total  Volume 

Load i ng 
Density 

W/V 

W/V 

1.20  psi 

0.7  25 
psi 

P 

Me t  hod 

(lbs) 

C  3 

(  lbs/ ft  ) 

1  3 

( lbs/ft  ) 

R  -ft 

R  -ft 

w 

ps  i 

l 

22040 

0.624 

0.211 

— o - 

2131 

— o - 

— 

2 

22040 

0.624 

0.21) 

— 

1  539 

— 

3 

22040 

0.624 

0.211 

1379 

— 

— 

4 

22040 

0.624 

0.211 

1092 

1586 

277 

5 

22040 

0.624 

0.21  1 

94  7 

1397 

322 

6 

22040 

0.624 

0.211 

1080 

1  569 

273 

1 

220400 

0.624 

0.499 

4590 

— 

— 

2 

220400 

0.624 

0.499 

— 

3952 

— 

3 

220400 

0.624 

0.499 

24  70 

— 

— 

4 

220400 

0.624 

0.499 

160  7 

2  3  34 

467 

5 

220400 

0.624 

0.499 

1356 

1999 

513 

6 

220400 

0.624 

0.499 

1617 

2348 

47  1 

NOTH:  Kxplosive  mass  and  chamber  volume  increased  10  times  and  total 

volume  increased  4  times. 
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3.2.3  Comparision  of  Six  Methods  -  Decrease  in  Tunnel  Diameter.  In 
table  4,  the  explosive  mass  was  increase*)  to  500000  lbs.  The  loading  density 
of  the  chamber  and  total  volume  remained  the  same.  The  only  difference  is  in 
the  diameter  of  the  exit  tunnel.  Here  you  can  see  that  Methods  1-3,  which  do 
not  use  the  tunnel  diameter  in  their  equations,  have  the  same  calculated 
distance,  while  Methods  4 -ft  show  a  reduction  in  distance  of  approximately  39%, 
which  corresponds  to  the  reduction  in  tunnel  diameter. 

Table  4.  Comparision  of  Six  Methods  -  Decrease  in  Tunnel  Diameter 


Charge 
Mass  W 

Loading 

W/V 

Densi ty 

W/V 

1 .20  psi 

0.725 

psi 

P 

Tunne  1 
Diameter 

Me  t  hod 

(  lbs) 

C  3 

( lbs/ f t  ) 

1  3 

( Ibs/ft  ) 

K  -ft 

— “O - 

K  -ft 
— o - 

w 

psi 

(ft) 

1 

500000 

6.24 

4.99 

6032 

— 

— 

16.6 

2 

500000 

6.24 

4.99 

— 

6790 

— 

3 

500000 

6.24 

4.99 

7609 

— 

— 

4 

500000 

6.24 

4.^9 

4519 

6561 

1888 

5 

500000 

6.24 

4.99 

3535 

5211 

1780 

6 

500000 

6.24 

4.99 

6  310 

9161 

2964 

1 

500000 

6.24 

4.99 

6032 

— 

— 

10.0 

2 

500000 

6.24 

4.99 

— 

6790 

— 

3 

500000 

6.24 

4.99 

7609 

— 

— 

4 

500000 

6.24 

4.99 

2722 

3952 

1888 

5 

500000 

6.24 

4.99 

2129 

3139 

1780 

6 

500000 

6.24 

4.99 

3800 

5518 

2964 
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3.3  Tunnel  Junctions.  When  there  are  two  exit  tunnels  nnH  they  are 
separated  enough  so  that  there  is  no  enhancement  between  them,  in  Equation  13 
the  value  of  n  becomes  2.  When  Method  3  is  used  to  calculate  values  in  Table 

3 

2  for  an  explosive  mass  of  11020  lbs  and  a  loading  density  of  0.312  Ibs/ft  , 

then  II  was  calculated  as  941  feet.  If  a  value  of  n  =  2  is  used  in  Equation 
o 

13,  the  distance  is  reduced  to  767  feet. 

It  is  suggested  by  the  author  that  a  new  approach  be  taken  when  there  are 
tunnel  branches  or  junctions.  This  new  method  would  reduce  the  transmitted 
pressure  by  factors  based  on  shock  tube  experiments.  These  reduction  factors 

are  presented  in  Figure  5  and  are  based  on  data  in  Reference  23.  The  90° 
tunnel  junction  data  do  not  follow  a  simple  equation  and,  therefore,  the  curve 
presented  in  Figure  f>  should  be  used. 

If  we  make  a  comparison  between  Methods  3  and  6  and  assume  a  Y  junction  in 
the  tunnel  system  that  gives  two  exit  tunnels  that  do  not  cause  any  exterior 
enhancement,  then  the  inhabited  building  distance  will  change  as  follows. 

Using  the  11020  lbs  in  Table  2,  the  distance  using  Method  3  is  941  feet.  If 
we  use  n  =  2  in  Equation  13,  this  distance  reduces  to  767  feet.  This  is  a 
reduction  of  approximate ly  Id?;. 

Now  using  Equation  7  to  calculate  I*  ,  we  find  equal  to  182  psi.  With  a 

Y  junction  as  shown  in  Figure  5c,  P  would  be  multiplied  by  0.65  to  become  UR 

w 

psi.  With  P^  equal  to  118  in  Equation  9,  the  inhabited  building  distance 
reduces  from  800  feet  down  to  581  feet.  This  is  a  reduction  of  27%. 

This  implies  that  using  n  *  2  may  be  conservative  and  that  the  941  feet 
should  reduce  to  687  feet  rather  than  767  feet. 

It  should  also  be  noted  that  with  a  tunnel  junction  as  shown  in  Figure 
5a,  there  would  be  different  exit  pressures  at  the  end  of  the  two  tunnels. 

The  inhabited  building  distance  would  also  be  different  in  front  of  the  two 
exits. 

The  reduction  in  pressure  propagating  through  the  different  junctions 
applies  only  if  the  tunnel  cross  sectioned  area  of  each  branch  remains  the 
same.  In  configurations  where  there  is  a  reduction  or  increase  in  the  cross 
section  area  of  the  tunnel,  then  these  conditions  should  be  treated  on  an 
individual  basis. 

24 

An  extensive  series  of  tests  were  conducted  by  Switzerland.  The  values 
given  in  Figure  5  compare  quite  well  with  the  results  reported  in  Reference 
24.  The  BRL  value  of  transmitted  pressure  of  0.80  P  in  Figure  5a  compares 

with  a  Reference  24  value  of  0.83  P  . 

s 
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90  Degree  Dead-end  Tunnel 


Y  Junction 
(C) 


Figure  5.  Transmitted  Pressure  versus  Input  Pressure 
for  Various  Tunnel  Junctions 


The  side  tunnel  (Figure  5n)  values  from  Reference  24  are  plotted  in  Figure 

6  along  with  the  BRL-generated  curve.  The  lM)°  dead-end  tunnel  (Figure  5b) 

value  for  the  transmitted  shock,  is  0.7  P  ,  while  the  Reference  24  value  shows 

s 

a  spread  of  0.57  to  0.68  P  ,  which  appears  to  be  partially  a  function  of 

incident  pressure.  The  Y  junction  transmitted  pressure  values  from  Reference 

24  for  equal  to  15°  through  90°  range  from  0.65  Pg  at  14.5  psi  down  to  0.58 

P  at  I'M)  psi.  This  compares  with  a  value  of  0.65  F  developed  in  Figure  5c 
s  s 

from  URL  data  in  Reference  23. 

4.  CONCLUSIONS 

4.1  Weaknesses.  It  is  impossible  to  establish  one  or  two  equations  that 
wild,  be  universally  accepted  and  that  fit  all  underground  storage  sites.  This 
report  has  presented,  discussed,  and  compared  the  results  of  six  methods 
proposed  for  determining  the  safe,  inhabited  building  distance.  All  methods 
have  certain  weaknesses,  some  more  than  others.  In  the  opinion  of  the  author, 
certain  parameters  should  be  known.  These  are  as  follows: 


Storage  Chamber  Volume 
Exit  Tunnel  Volume 
Loading  Density 
Explosive  Distribution  and 
Containment 

Chamber  Pressure  for  Specific 
Explosives 


Chamber  Diameter 

Tunnel  Diameter 

Tunnel  Junctions  (If  any.) 

Tunnel  Roughness 

Terrain  Outside  of  Tunnel 


All  of  these  variables  will  affect  in  some  way  the  overpressure  propagated 
outside  of  the  tunnel.  One  other  variable  not  dealt  with  is  the  location, 
confinement,  and  point  of  initiation  of  the  explosive  source.  The  major 
portion  of  scaled  model  tests  has  been  conducted  with  linear  charges  placed 
along  the  centerline  of  the  chamber  or  near  spherical  charges  placed  near  the 
entrance  to  the  storage  chamber.  When  in  a  real  storage  scenario,  there  will 
be  pallets  and  boxes  of  munitions  stored  throughout  the  chamber  and  on  the 
floor.  Most  of  the  munitions  will  have  some  kind  of  containment,  from  the 
thin  skin  of  rocket  motors  to  the  thick  casing  of  general  purpose  bombs.  The 
effect  of  containment  on  the  build-up  of  gas  pressure  within  the  storage 
chamber  has  not  been  fully  addressed. 

4.2  Recommendations.  It  has  been  shown  that  Methods  4  or  6  give  the  most 
consistent  values,  and  the  inhabited  building  distances  vary  only  a  few 

percent  in  the  medium  loading  densities,  i.e.,  less  than  0.624  Ib/ft  .  At  the 
higher  loading  densities,  it  is  recommended  that  Method  b  be  used  in  any 
prediction  calculation.  It  can  he  seen  in  Figure  2  that  using  the  loading 
density  (W/V^)  as  an  input  parameter  in  Equation  4  will  give  different  exit 

pressures,  than  using  the  static  pressure  (P  ),  which  is  based  on  (W/Vt)  as 

the  input  parameter  in  F.quation  7. 
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Having  available  this  list  of  ten  variables,  thete  is  still  no  assurance 
that  a  precise  prediction  can  be  made.  The  methods  presented  here  should  be 
used  as  guides  and  not  for  planning  and  construction  of  new  sites.  When 
planning  the  location  of  a  new  site,  it  is  recommended  that  a  scaled  model  of 
the  site  be  constructed  and  tests  conducted  to  determine  the  range  for 
inhabited  buildings.  This  is  also  true  where  there  may  be  a  controversy  over 
a  specific,  existing  site. 
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